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Abstract. The role of reflection symmetry breaking for the character of the appearance of replica symmetric 
spin glass state is investigated. We establish the following symmetry rule for classical systems with one 
order parameter in the replica symmetric mean field approximation. If in the pure system the transition 
to the ordered phase is of the second order, then in the corresponding random system the glass regime 
appears as a result of a phase transition; if the transition in the pure system is of the first order, then glass 
and ordered regimes grow continuously in the random system on cooling. 
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The crucial role of the reflection symmetry for the 
character of phase transition in nonrandom mean-field 
(MF) models is well known (see, e.g., the textbook Q). 
Generally speaking the presence of the terms without re- 
flection symmetry causes the first order phase transition, 
while in the absence of such terms the transition is of the 
second order. Usually this result is obtained in the frame 
of the phenomenological approach based on the Ginsburg- 
Landau (GL) effective Hamiltonian that can be easily ob- 
tained for any Hamiltonian through the Hubbard-Stratanovich 
identity for the partition function. On a slightly "more 
microscopic level" the statement about the order of phase 
transition can be demonstrated considering the MF one- 
site equations as is done below (for continuous case see 
the Eqs. (11-13) of §). 

Let us consider a system of particles on lattice sites i, j 
with Hamiltonian 



fluctuations ip. Using the Hubbard-Stratanovich identity 
the partition sum can be written in the following form: 



Z = Tr exp 
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2^ 



det (pjy 
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[2tt 
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J[d<p t exp{-f3H eff } (3) 



H = - - ^2 JijUiUj, 



(1) 



where U is a diagonal operator with TrJ7 = 0, the inter- 
actions Jtj being such that the MF approximation gives 
exact solution. For example, Jy may be random exchange 
interactions with Gaussian probability distribution 



P{Jij) = 



1 



: exp 



(Jy — Jo) 



2J 2 



(2) 



with J = J/yN , Jo = Jo/N. In nonrandom pure case 
J = 0. Let us construct GL Hamiltonian as an integral on 
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The effective Hamiltonian contains the powers n of the 
fluctuations ip for which Tr(U n ) is nonzero. 

When dealing with random MF systems the role of cu- 
bic term is discussed usually in connection with the replica 
symmetry breaking (RSB) solution of the equations for 
glass order parameters. The absence of reflection symme- 
try causing cubic terms in GL Hamiltonian for regular 
nonrandom system results in a special form of RSB free 
energy functional for random case, too. This form gives 
the discontinuity for the order parameter, the stability of 
the first stage RSB solution and other specific features 
(see, for example, Refs. J|,[|,[|). However, in this paper 
we accent our attention on the replica symmetric (RS) 
approach. We investigate the role of the reflection symme- 
try for the behavior of RS solution for spin-glass-like MF 
systems with one regular order parameter in the Hamilto- 
nian and formulate a kind of symmetry rule for the type 
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of the growing of glass regime: if in the nonrandom (pure) 
system the transition to the ordered phase is of the sec- 
ond order, then in the corresponding random system the 
glass regime appears as a result of a phase transition; if 
the transition in the pure system is of the first order, then 
in the random system the glass regime grows continuously 
on cooling. In fact, one can imagine that both the first or- 
der phase transition in nonrandom systems as well as the 
continuous growing of the glass regime in random systems 
are caused by some kind of internal fields appearing due to 
the algebra of operators U . Some previous results can be 
found in Q. It is well known that a number of magnetic 
compounds with phase transition to spin-glass state be- 
haves qualitatively as the Sherrington-Kirkpatrick (SK) 
model JtJ (without ip 3 term). We would like to bring the 
readers attention to the fact that the case with the con- 
tinuous growing of the glass regime (with ip 3 term) also is 
a reality and exists in some mixed crystals 

Now let us consider the system (|l|)- (^|). Using replica 
approach (see, e.g. [Q) we can write the free energy in 
the form 



(F) j/NkT = - lim - max! - V 
n— >o n y z — ' 



(x a ) 2 



lnTr|;yc.} exp 




J2 y ap tU a U 

a>(3 
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where 




q a0 = (y «/3 r tr A = {U <*U0) eS . 



(5) 



as follows from the saddle point equations. Here t = J/kT 
and averaging is performed with the effective Hamiltonian 



-n eff = ]T ^™ a u a +Y, \p a {u a ) 2 +Y. t 2 ^ p u a uP 

a a ct>(3 



In the RS approximation the free energy has the form: 



F = -NkTi- ( — ) — +t 2 — -t 2 —- 
] \ kT 2 4 4 



f 



C L= exp | I In Tr 
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(e)}} (6) 



Here 



ztyfq + m 




u + t 2V -^-u 2 . 



The order parameters are: m is the regular order param- 
eter (an analog of magnetic moment in spin glasses), q is 
the glass order parameter and p is an auxiliary order pa- 
rameter. The extremum conditions for the free energy (|]) 
give the following equations for these order parameters: 
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In general case the order parameter p is not indepen- 
dent from the others. For Ising spin glass § U = S, 
(U) 2 = (S) 2 = 1 = p and Eq. (^) reduces to an iden- 
tity. For the quadrupolar glass with J = 1 KM we have 

~ 2 

U = Q = 3J Z — 2, so that p — 2 — m. However, for the 
quadrupolar glass with J = 2 JlTJ] p is independent order 
parameter: the subalgebra contains the operators Q, (Q) 2 
and the unit matrix E, so that only (Q) 3 is not indepen- 
dent and can be represented as a linear combination of Q, 
(Q) 2 and E. In the case of S = 1 spin glass JllJ] we have 
S z = 0, ±1, p is independent parameter and has the physi- 
cal meaning defined by the equality S 2 — 1/3(2 + Q). Now 
p = (1/3) (2 + 7712), where m 2 is the quadrupolar regular 
order parameter. 

In the random case Jo = the high temperature ex- 
pansion of the equations (|?|) - (|J) has the form: 

t 2 t 4 t 4 

m = —A 3 p + —p 2 {A 5 - 2A 2 A 3 ) - —A 2 A 3 q 2 - 



t 2 A\q- 



--A 2 A 3 pq + ... (10) 
-A\p 2 + %q 2 {A\ - AA 3 ) + t 4 A 2 A A qp + ... 



(11) 



P = A 2 + jp(A± 
t 4 



A l) + zjP 2 (A - ZA 2 A 4 + 2A 3 )- 



t 4 



-- q 'A 2 (A-4-Ai)--pqAi + .. 

If m = 0, q = and p / then 

Tr[f7exp(0)] = Tr[U exp(i 2 p[/ 2 /2)] 



(12) 



(13) 



is a sum of A n = Tr(C/™) with odd n only. So if A^+i) = 
then Eqs.(f7|) - (0) have the solution m = 0, q = 0. If 
^4(2fc+i) 7^ Eqs.(|7]) - (^) have no trivial solution m = 0, 
q = at any temperature. The high temperature expan- 



sion gives explicitly m ~ t A 3 ,q 
parameter grows continuously on cooling. 



t 4 A 2 . 



The glass order 



N.V.Gribova, V.N.Ryzhov, T.I.Schelkacheva,, E.E.Tareyeva: Reflection symmetry in mean-field replica-symmetric 



3 



In the case A^k+i) = one can obtain from (||) the 
bifurcation point T* where the nontrivial solution for the 
glass order parameter q appears under the condition m = 0. 
The equations for T b have the form: 

p(t c ) = l/t c ; t c = JjkT h c 

p(t c ) = Ti[U 2 exp{t 2 p(t c )U 2 /2)}. 

The replica symmetric solution is stable unless the 
replicon mode energy A is nonzero. For our model we have: 
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The character of the replica symmetry breaking is defined 
by the term az{8q a pf in the RSB free energy (Q) expan- 
sion near RS free energy (^|). We have 
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then Tr 
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and g = at 



T = T c so that 03 — and at T c the full replica symmetry 
breaking according to the Parisi scheme takes place. If 
A-2n+i 7^ then a3 ^ and the first stage RSB solution 
is stable. 

In the regular case (J = 0, Jo ^ 0) we have from 
Eqs.(|) - (|) q = m 2 and 
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= A, = Ao, A = A + r. 
0, the solution of ( |l7j ) has the form 



A 2 



%tA\ 



^A 2 Ao z 



iA 4 A 3 



3A 2 



A 4 



(18) 



and we have the second order phase transition. If A3 7^ 
the solution of Eq. (n5f) in the vicinity of T b is 



2rA| 
A3 



(19) 



In this case the phase transition is of the first order. The 
transition temperature can be obtained from the compari- 
son of free energies of ordered and disordered phases, (see, 
e.g. Ref-ill). 

So one can establish the following symmetry rule for 
classical systems with one (pure) order parameter in the 
replica-symmetric MFA. If in the pure system the transi- 
tion to the ordered phase is of the second order, then in the 
corresponding random system the glass regime appears as 
a result of a phase transition; if the transition in the pure 
system is of the first order, then in the random system 
glass and ordered regimes grow smoothly on cooling. 

As an example let us now consider the Hamiltonian 

(@) with U = Q + r]V, Q = 3J Z 2 -2, V = V3{J X 2 - 
-2 

J y ), J = 1, J z = 1,0,-1. This model describes random 
quadrupole interaction. A molecular quadrupolar moment 
is the second-rank tensorial operator with five independent 
components. In the principal axes frame only two of them 
remain: Q and V. It is easy to show that Q 2 = 2 — Q, 
V 2 = 2 + Q, QV = VQ = V. Here 77 > is a tuning 
parameter. 



(Q + vV) 2 = 2(1 + if) + 2 V V + (1 - r] 2 )Q. 



(20) 



The corresponding Ginzburg-Landay Hamiltonian writ- 
ten in terms of fluctuation fields (pi has the form: 

P H *ff = ^fiW^Vi — Nln3— (21) 

ij 

[pta + V 2 ) + tfiv 2 1/3) + <pf(l + ^ 2 ) 2 /4 + . 

i 

The RS free energy has the form 



N 



These mean field equations have the trivial solution m = 
andpo — Tr(£/ 2 )/Tr(_E) for all temperatures. A non-trivial 
solution for m appears at the bifurcation point kT c b / Jq = 
p Q . Near T c b the expansion of (|l5|) - (|l^) gives the bifur- 
cation equation 



F 



-NkT\- 

I \kT 2 4 



_ t 2^2 +t 2(l + r] 2 ) 2 -t 2 (l + r 1 2 )q + + 



dz 



exp 



blip} (22) 



Here 



ip = exp(-2i?i) + exp(z?i) [exp(tf 2 ) + exp(-i? 2 )] 
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0! = *Vff + mi (^) + *V"1) 



j? 2 - r?V / 3t 2 [2(l+r; 2 )-(m 2 + g )]+77V3 



W ; 2 2 



m i 



fcT 



The order parameters are m\ Q + r\V ^$>, q is the 
corresponding glass order parameter, the order parameter 
m 2 ^<C (1 — y 2 )Q — 2r\V ^> appears due to the operators 
algebra (see Eq. (ppf)) in close analogy with the appearing 
of quadrupole order parameter in the case of spin glass 
with S = 1. We use m 2 instead of p: 

p = 2(if + 1) - m 2 . 

The equations for the order parameters have the form: 



mi = 



dz 
n/2tF 



exp 



m 2 = 



dz 
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ex P 7T 
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«,9, t '/V o fl , 
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(25) 



It is easy to see that 

(Q + r/F) 3 = 2(3t7 2 - 1) + 3(1 + V 2 )(Q + r,V)}; 

Tr(Q + rjVf = 6(3r? 2 - 1) 

Using these equalities and Eqs. (|2l|) one can show that 
the behavior of the system under consideration is quite 
different in the cases rj = 1/V3 and rj ^ 1/V3. 

Let us consider first the re gula r case (J = 0, Jo 7^ 0). 
Now we have from Eqs. (|l5|)- fll^): 



mi 



16 /" Q 2 



m 1 



*r(l 



3t/ 2 



15*0. 



The case 77 = 1/V3 has transparent physical meaning for 
Jo > 0, second order transition for spin one ferromagnet 
being an example. In the case 3^ 2 - 1 / we obtain the 
first order phase transition for 7/ > l/v3 and Jo > or 
for 77 < l/v3 and Jq < 0. It should be noted that an 
example of such a behavior is presented by the case 77 = 
describing the orientational quadrupolar ordering in o—H 2 
or p — D 2 (see ) . We can re write equation ( |l5| ) as 

A) = (y - 4A) exp(3y/2) + 2y + 4A = (26) 



where y = -2m (-§), A = (#,). 

The bifurcation points A^ for <£(y, A) |lJ 
from the condition I J =0. We have 



are obtained 



A x = l/2;pi = 0;m ( l)=0; 

A 2 = 1/2.18; yi = 0.7; m { 2) = -0.8. 

The form of the solutions can be obtained from the 
expansion of (^6|) in powers of r = A — \ a and £ — y — y a . 
Because TrQ 3 ^ in the vicinity of T b = J //cAi we have 
£ ~ r. However in the vicinity of T 2 b = 1.09Ti b we have 
£ ~ V^- The physical solution near T 2 b takes place only 
for /1 > 0, T < T 2 b . It is a turning point. The sharp phase 
transition with the jump of m is between T\ and T 2 . 
The transition point is calculated from the comparison of 
energies. Other cases with 7/ 7^ l/\/3 can be described in 
a similar way. 

Let us consider now the random case (J 7^ 0,Jo = 
0). The first terms of the high temperature expansion of 
Eqs. (@) - (||) in the case U = Q + r/V can be easily 
obtained from (noh - (|12|) and have the form: 



mi = 2t 2 ab + /j 4 a6 3 — t 2 am 2 — t A m 2 ab 2 H m 2 ab— 

4 

-U A ab 2 q - 2t 4 afeg 2 + 2t i abm 2 q, (27) 
g = 4tVg + 4< 4 a 2 6 2 - 4t 4 6a 2 m 2 
+/j 4 a 2 m 2 + 8£ 4 o 4 g - U i b 3 qm 2 + 2t i q 2 (a 2 - 8b 2 ), (28) 
m 2 = -2bH 2 + t 2 b 2 m 2 - t% 2 {2a 2 - b 2 ) + t 4 bm 2 (2a 2 - b 2 ) 
t 4 

- — (2a 2 - b 2 )m\ + 2i 4 6 3 g 2 + 4t 4 a 2 bq - 2rj 4 a 2 m 2 g. (29) 



Here 



1. 



l + T) . 



If r\ = l/-\/3 the Hamiltonian becomes analogous to 
the S = 1 spin glass Hamiltonian. There is the trivial so- 
lution q = when T > T c b because A3 = 0. The solution 
q 7^ appears at the point kT c b /J rss 3.2 (for m = 0) |Q . 
The specific heat function loses smoothness at this point 
(see fig.l). In the case a 7^ there is no trivial solution of 
Eq. (28) and g, m 2 and mi grows continuously on cool- 
ing. This case is analogous to quadrupolar glass with two 
internal fields. The absence of the zero solution at high 
temperatures for glass order parameter is unambiguously 
connected with the absence of reflection symmetry. 

In the figures we present the behavior of the order 
parameters for different cases as well as the specific heat 



Cy 




kN d(kT/J) 



Jo \ dmi 

— mi — 

J / d{kT/J) 



(30) 



Let us note that the symmetry rule for the type of 
the growing of glass regime and phase transition in pure 
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Fig. 1. Order parameters for r\ = 0.5/\/3 1) Jo = 0; 2) 
Jo/ J = 1.5; and 3) specific heat for a) Jo/ J = 0; b)Jo/J = 1.5; 
c)J /J = 2. 
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Fig. 2. Order parameters for r\ — l/v3 1) Jo = 0; 2) Jo/ J = 
1.5; and 3) specific heat for a) Jo/ J = ; b)Jo/J = 1.5. 



systems have been obtained in MFA. It is interesting to 
notice that in the Bethe approximation (and in some other 
cluster approximations) the internal fields appear in the 
case of nonzero cubic terms and cause a smooth changing 
of the order parameter without phase transition in non- 
random system if in MFA the phase transition is of the 
first order (see, e.g., In the systems with reflection 

symmetry and the MFA second order phase transition the 
rather sharp phase transition from nonzero order param- 
eter phase remains in the Bethe approximation. In some 
way the replica approach to random systems effects anal- 
ogously to Bethe approximation. 



The simpliest way to demonstrate this fact is to con- 
sider cluster approximation of |[6| |. The main equation 
is the selfconsistency condition given in terms of effective 
mean field ip produced by the neighbors and acting on the 
particle on each site of the cluster. 

For the Hamiltonian ([j]) in nonrandom case with near- 
est neighbor interaction we obtain the following equations 
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<u>= 



Trlqexp -/lffi(V>) 



Trexp 



TrJ7i exp 



Trexp -0H a y>) 



(31) 



where s is the number of particles in the cluster. For the 
simplest case s = 2 we have where Hi = — U\tp\B.2 = 
- JUiU 2 - ip (1 - I/7) (t7i + U2) ; is the operator U 
on i-th site; 7 is the number of nearest neighbors. The 
effective mean field ip can be obtained from (pl|). The 
high temperature expansion for ( ^l| ) is 



V>/3JTr (c/ 2 ) + ^ 2 ^-Tr (l> 3 ) = 7/;a/3Tr (l> 2 ) + 



/3 2 J 2 



■Tr 



(l> 2 ) Tr (l> 3 ) + ^ 



2 02 



■Tr 



Tr 



(32) 



where cv = 1 — I/7. 

It is easy to see from ( p2[ ) that there is no solution 

ip = for the case Tr *j 7^ 0. So there is no phase 
transition and the order parameter grows smoothly on 
cooling. When Tr \ = a bifurcation point can ex- 
ist for Eq. ( |32| ) and in the pure system we obtain a phase 
transition. 

In conclusion we formulate a kind of symmetry rule 
for the way of the appearance of glass regime: if in the 
nonrandom (pure) system the transition to the ordered 
phase is of the second order, then in the corresponding 
random system the glass regime appears as a result of a 
phase transition; if the transition in the pure system is of 
the first order, then in the random system the glass regime 
grows continuously on cooling || . 
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